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ABSTRACT 

Recent high-resolution and high-cadence observations have surprisingly suggested that prominence 
barbs exhibit apparent rotating motions suggestive of a tornado-like structure. Additional evidence 
has been provided by Doppler measurements. The observations reveal opposite velocities for both hot 
and cool plasma on the two sides of a prominence barb. This motion is persistent for several hours 
and has been interpreted in terms of rotational motion of prominence feet. Several authors suggest 
that such barb motions are rotating helical structures around a vertical axis similar to tornadoes 
on Earth. One of the difficulties of such a proposal is how to support cool prominence plasma in 
almost-vertical structures against gravity. In this work we model analytically a tornado-like structure 
and try to determine possible mechanisms to support the prominence plasma. We have found that 
the Lorentz force can indeed support the barb plasma provided the magnetic structure is sufficiently 
twisted and/or significant poloidal flows are present. 


1. INTRODUCTION 

Recent high-resolution and high-cadence observations 
have revealed a possible rotating motion of prominence 
barbs around a vertical axi s. These have been called 
barb tornadoes (lPriestll20I4fl due to the similarity of the 
apparent motion of such stru c tures on the limb to ter¬ 
restrial tornadoes. iSu et al.l (j20I2t) reported an event 
where a prominence shows an apparent rotating motion 
with velocities of up to 8kms“^. The observation shows 
cool prominence plasma seen in absorption in the I7I 
A passband of the SDO/AIA instrument. The authors 
argue that the motion is a rotation projected on the so¬ 
lar limb plane. However, such projected motions are 
also compatible with osc illations and counter-st reaming 
flows, as pointed out bv iPanasenco et iaLl (j2014l l. More 
events, up to 201 barb-t ornadoes, have been reported by 
iWedemever et al.l (|2013ll using AIA 17lA images. How¬ 
ever, these authors studied mainly the morphology of the 
b arbs and deduced a wide ra nge of sizes and lifetimes. 

lOrozco Suarez etldl (|2012h reported Doppler shifts us¬ 
ing He 11083.0 nm triplet data from the German Vacuum 
Tower Telescope (VTT) of the Observatorio del Teide 
(Tenerife, Spain). The observations revealed opposite ve¬ 
locities at the edges of the prominence feet of ± 6 kms“^ 
along a slit placed close to the solar surface suggesting ro¬ 
tation of the cool prominence feet. The width of the feet 
is about 20" indi cating a period of rotation of 4 hours. 

More recently, iSu et al.l (j2014ll also reported Doppler 
shifts in a prominence pillar using the Fe XII 195A line 
of the EIS instrument onboard the Hinode satellite. The 
observations revealed a bipolar velocity pattern along 
the whole vertical prominence pillar. The velocity is al¬ 
most zero at the tornado axis and increases linearly up to 
^ ±5kms“^ at the two edges of the observed structure. 
This indicates that also the million-degree plasma related 
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to the tornado-like prominence may be rotating. The 
EUV bands of SDO/AIA reveal that the cool plasma seen 
in absorpti on moves in consort with t he hot plasma. Ad¬ 
ditionally, iMartmez-Gonzalez et al.l (|2015f) . have found 
evidence of helical magnetic structure simultaneously at 
two prominence feet. All this evidence suggests that barb 
tornadoes are rotating vertical structures, nevertheless 
more observational evidence is needed. 

The existence of such structures opens new questions 
concerning the origin of the tornado rotation and the in- 
fluence of the rotating barb on the rest of the filam ent. 
IWedcmever-Bohm et al.l (|20I2f) and ISu et al.l (|20I2ll pro¬ 
posed that the barb-tornadoes are driv en by photo- 
spheri c vort ex flows of the kind observed bv iBrandt et al.l 
(1198811 and IBonet et al.l (1200811 : according to this view, 
the barb field lines are rooted in the vortices and the lat¬ 
ter’s rotating motion is transferred to the barb. The 
authors also proposed that the barb tornadoes inject 
chromospheric plasma and helicity into the upper fila¬ 
ment throughout the rotating barb. However, such vor¬ 
tex flows have yet not been observed at all below barbs, 
let alone as a matter of course. 

Another possibility for the origin of spiral motions is 
three-dimensional reconnection (at or above the pho¬ 
tosphere) during cancellation of photospheric magnetic 
fragments, since such re connection w ill naturally produce 
vortex motions (e. g.. IPriestI 1201411 and could also fuel 
a prominence with mass. A third possibility arises from 
the fact that a prominence represents a concentration of 
magnetic helicity in twisted magnetic fields. Thus, if part 
of a prominence dips down towards the photosphere, it 
is possible that such magnetic helicity and its associated 
twisting motions may be focused in the dip. 

In this letter we discuss how the massive cool plasma 
is supported against gravity in a helical magnetic struc¬ 
ture. We find that the barb-tornadoes bear many sim¬ 
ilarities to astrophysical plasma jets in which magneto¬ 
centrifugal forces accelerate the plasma. By using recent 
current tornado data and typical barb parameters we 
conclude that it is actually possible for the magnetic force 
to support and accelerate the cool barb plasma against 
gravity provided the structure is highly twisted and/or 
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significant poloidal flows are present. 

2. THE BARB-TORNADO MODEL 

Inspired by the observations, we consider in the fol¬ 
lowing an axisyminetric model rotating around a verti¬ 
cal axis. We use cylindrical coordinates (r, (f, z) with 
z coinciding with the rotation axis. The axisymmetry 
condition implies that all quantities are yj-independent: 
dj dip = 0. An important restriction on the axisymmet- 
ric magnetic field is that the field lines should become 
vertical as we approach the rotation axis (r —>■ 0). The 
observed structure appears to be stationary with no im¬ 
portant changes of shape in the EUV images and Doppler 
pattern for several hours. The Alfven speed is around one 
thousand kms“^ in the corona and of order one hundred 
kms“^ in the cool prominence plasma. Thus, the travel 
time of a magnetic perturbation is less than a minute 
along the vertical axis of the tornado. This indicates 
that, during the few hours the observation, the system 
has plenty of time to relax and produce a stationary 
magnetic structure. We therefore set djdt = 0 in the 
equations. In this situation the MHD induction and mo¬ 
mentum equations become, respectively, 

0 = V X (17 X B) , (1) 

0 = —p(i7-V)?/—Vp-I-—(V X B) X B-I- /95 , (2) 

Mo 



Fig. 1.— Schematic picture of a torna do-like mag n etic s tructure 
using the elementary force-free field of ISchatzmaru 1I1965IV Solid 
lines are the three dimensional representation of the magnetic field 
lines. Dashed lines are the poloidal field lines in the plane = 
constant. 
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where p, 17, B, p and g = —gCz are the plasma density, 
the velocity, magnetic field, gas pressure and gravity re¬ 
spectively. Both the plasma velocity and magnetic field 
can be naturally decomposed into poloidal and toroidal 
components, 

v = vp + v^e^, (3) 

B = Bp + B^ , (4) 

with Vp ■ = Bp ■ = 0. Only for illustration purposes 

we are showing i n Figure [1] the force-free solutions of 
ISchatzmai] (|1965ll . The poloidal planes (r, z) correspond 
to ip = constant surfaces. Given the axisymmetry, we 
can easily define the angular velocity D(r, z) using = 
r f2(r, z) e^. The advection (or inertial) term (y ■ V)l7 can 
be written as 

iv ■ V)v = (vp ■ V)z7p — rfl^ Or + — ■ V(r^D)e„ , (5) 

P 

i.e., advection in the poloidal plane plus two accelera¬ 
tion terms associated with the prescribed rotation profile 
0(r, z), the first of which is clearly the centripetal force 
of a circular motion. 

In a stationary regime there must be force balance in¬ 
cluding the inertial terms. For ease of notation, we use 
the symbol Fnm for the sum of pressure gradient, gravity 
and centrifugal force: 

Fnm - Vp - pge^ P fyrV?ep . (6) 


+ 


Mo 


f^pol 7 


( 8 ) 


where e_|_ are the unit vectors in the poloidal plane 

parallel and perpendicular to Bp, respectively, and Kpoi 
and s are the curvature and arc-length along the poloidal 
field lines (the latter is chosen such that Cy points in the 

sense of growing s illustrated as dashed lines in Fig. [T]). 
For future reference, we also note that Vp and Bp are 

defined through Vp ■ ey and Bp ■ ey, respectively, i.e., they 
can be positive or negative. Equation © represents the 
steady flow of plasma along field lines and provides clues 
concerning the support of the cool barb plasma; Equation 
([8]) represents the global transverse equilibrium of the 
magnetic structure. The remaining component of the 
equation of motion, the azimuthal component, is 

p(Fp-V)(r2D) = (ri?^). (9) 

As an additional remark concerning Figure [U note 
that in a barb-tornado structure the density along the 
poloidal field lines should be larger close to the axis than 
in the oute r regio ns of the structure. The observation 
of ISu et ahl (120141) . which shows a dense column in ab¬ 
sorption tapering off with height, may well be compatible 
with a helically opening structure as shown in the figure. 


The poloidal part of the equation of motion ([2) can then 
be split into into its components parallel and perpendic¬ 
ular to Bp, 


p dvl ^ p ^ _1 d{rB,p)^ 

2 ds 2/ror2 ds 


(7) 


2.1. Pure Rotation 

In a situation of pure rotation there is no plasma flow 
along the poloidal field {vp = 0), and the inertial term 
([5]) reduces to the centripetal component —prft'^e^. In 
this case, from Equation ([9]), we see that the azimuthal 
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component of the Lorentz force is zero, i.e., rB^p is con¬ 
stant along each poloidal field line. On the other hand, 
using the induction equation one can easily see that in 
this case id is constant along each give n field line, i n 
agreement with Ferraro’s isorotation law (lFerrarolll937f) . 
The magnitudes of rB^, and id are then determined by 
the boundary conditions of the problem. Given the con¬ 
stancy of rBtp, we see from the longitudinal equilibrium 
equation 0 that the Lorentz force has no longitudinal 
component, and so the equation reduces to 

Km • e|| = 0 . (10) 

The Lorentz force is thus purely poloidal and perpen¬ 
dicular to Bp. Calling 0 the angle between the poloidal 
field and the horizontal direction, equation (HOD may be 
written 

r fd^ COS0 — - ^ — g sin0 = 0 , (11) 

p as 

implying a balance of the non-magnetic forces. Note, in 
particular, that the centrifugal term could help support 
the plasma against gravity if the field lines are sufficiently 
close to horizontal (i.e., sin6( > 0 sufficiently small). 

Is the foregoing purely-rotating stationary equilibrium 
a realistic possibility for the obse rved apparent barb tor¬ 
nadoes? The latest observations (|Su et al.ll2014ll reveal a 
rotational velocity of 5 km at the edges of the struc- 
ture with r = 2" = 1.5 Mm, so Id 3 x 10“^ rad s“^. 
iWedemever et al.1 l)2013ll found a verage val ues of barb 
widths of r Ri 2" in agreement with lSu et al.L With these 
values we can compare the centrifugal acceleration to the 
gravitational acceleration, 

r 

-« 0.06 = tan(3.4°) . (12) 

9 

With the values found by lOrozco Suarez et al.l (|2012[) 
this ratio is even smaller. The only way to have centrifu¬ 
gal support of the plasma is then for the field lines to 
be almost horizontal, which contradicts the tornado pic¬ 
ture. Can a pressure gradient help support the plasma 
in a non-horizontal field? Hot, coronal plasma, can be 
supported by a pressure gradient against gravity across 
coronal distances even in vertical field lines. However, 
for cool prominence plasma the pressure scale-height is 
too small to balance gravity in barbs as tall as those ob¬ 
served. Even if the field lines close to the rotation axis 
were filled with hot coronal plasma, the magnetic struc¬ 
ture would have to turn almost horizontal (say, 3°, as 
given in Eq. HI in the region holding cool prominence 
plasma. We conclude that the centrifugal force asso¬ 
ciated with the rotation cannot support the cool barb 
plasma against gravity in a helical field structure. 

Another way to illustrate this conclusion is to estimate 
the rotational speed necess ary to ha ve purely centrifu¬ 
gal support, namely v^p = ^rg tan 9 . Assuming a mag¬ 
netic field inclination of 45°, say, at the edges of the ob¬ 
served barb tornadoes, the rotation vel ocity wou l d have 
to be 20 km s“^ if we use the data of ISu et HI (|2014f) 
and 4 5 km s“^ when using those of lOrozco Suarez et al.l 
(|2012f) . much larger than the measured speeds. 

2.2. General case with poloidal flows 


We consider now a more general scenario allowing for 
flows in the poloidal direction. Now the simple situ¬ 
ation of constant H and rBp along the field lines no 
longer applies. Checking for instance the azimuthal 
equation of motion dH), we see that a change of spe¬ 
cific axial angular momentum of the plasma elements 
associated with the poloidal motion must be associated 
with a non-zero toroidal component of the Lorentz force, 
{Bp/r) ■ \7{r Bp) Bp,. Hence, in general rBp can no 
longer be constant along field lines. In this situation, 
the Lorentz force also has a non-zero projection along 
poloidal field lines (Eq. [T]), 


def 1 d{rBp)'^ 
Mlong--2^^p2 Qg 


(13) 


In spite of the added complication of this general sce¬ 
nario, there is a set of conserved quantities along the field 
lines. One can obtain them by following the gene ral pro¬ 
cedure used in the theory of astrophysical jets (IMestell 
119611 : iLovelace et al.l 119861) . The induction equation re¬ 
quires that V X (Fp X Bp) = 0, so the poloidal velocity 
must be parallel to the poloidal field lines, for otherwise 
a singularity would arise in the toroidal electric field at 
the 2 ;-axis. So, we can write 


Vp = K{r, z) Bp . (14) 

Integrating the MHD equations along field lines and 
simplifying the resulting expressions, a set of conserved 
quantities results, namely. 


PoP^ = > 

(15) 

= w , 

PoPr 

(16) 

K 

(17) 


where K, W, and A are all constant along each poloidal 
field line. Those relations allow us to find an explicit 
expression for rBp along the poloidal field. One can write 
it in terms of the poloidal Alfven Mach number 


M\ 




Bj/PoP PoP ’ 
and of the Alfven radius^ ta^ defined by 

9 def A. 

as follows: 


rB^=KW 


I-Ml 


(18) 


(19) 


( 20 ) 


Expression (j20|) become singular when Ml^ —1. In the 
classical wind solutions, the flow speed increases from 
sub-Alfvenic to super-Alfvenic at a given point, and, to 
avoid the singularity, this transition must happen pre¬ 
cisely at the Alfven radius, r = r a- In such models, this 
requirement serves as an internal boundary condition to 
pick up the desired trans-Alfvenic solution instead of the 
a completely sub-Alfvenic or completely super-Alfvenic 
solution. In our case, however, it is unlikely that the 
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flows reach Alfvenic values within the tornado; rather, as 
shown below, we expect M\ to remain below 1 through¬ 
out. 

Following Eq. ([7]) and the results of Section lTAl to sup¬ 
port the cool plasma along field lines the magnetic force 
Fmiouq of Ed- should be larger than the centrifugal 

force, fc = r. From Eqs. (USD and (El), the ratio of 
these forces is 


F, 


Mlong 


fc 



( 21 ) 


in terms of the ratios (Tb and T„) of azimuthal to 
poloidal components for the magnetic fields and plasma 
velocities: 


d^B^p df^Vp 

1 Q - Sriicl 1 y - , 

Br, Vr, 


and the length Lq, defined as 


Lq = 


91n(r^f7) 


ds 


-1 


( 22 ) 


(23) 


The characteristic length (1751) can be assumed to be of or¬ 
der the radial size of the system, i.e., Lq ~ 0(r). Hence, 
from Equation ED, it is necessary for the value of the 
parameter \TB/Ty\ to be large in order to have a barb- 
tornado dominated by the magnetic force. We see that 
a high level of magnetic twist, and a comparatively large 
poloidal velocity (compared to the rotational velocity) 
are necessary to provide magnetic support. 

We can derive a more exact estimate of the possibilities 
of magnetic support of tornado plasma against gravity 
using expression (113)) for the longitudinal magnetic force 
combined with expression (l20)) . First of all, we can use 
the definitions of W and K to write their product in the 
form 

KW = Ml ^{Ty-TB) (24) 

p p 


keeping in mind, nevertheless, that A", W (and r^) 
are constant along each field line. We can now derive 
Eq. (l20l) with respect to the arc-length and, after a little 
algebra, obtain the general expression 


F, 


M long 


2 Bp Bp 
r po 1- 


cos 9 


(i + f)T,-r,] 


(25) 

where we have used cos0 = dr/ds, and called a the scale 
of variation of p along poloidal field lines in terms of the 
cylindrical radius, as follows: 


de/ r 9 log p 
cos 9 ds 


(26) 


For the estimate that follows, it is best to write Equa¬ 
tion (1251) in terms of an observed quantity, namely Vp, 
and of the ratio TB/Ty, which turns out to be the essen¬ 
tial dimensionless variable in the resulting expression. 
Further, to calibrate the possibilities of magnetic sup¬ 
port, we normalize (1251) with respect to gravity along the 
poloidal field, pg sin0. The result is: 


Em long 

pg sin0 


2 



Tb\ 

TyJ 




Ie 

Ty 


- 1 


cot6». (27) 


For the applications to prominence barbs we envisage 
cases with < 1. The term 1 — in Eq. [77] 
therefore favors magnetic support. In fact, in many cases 
<C 1 and so for simplicity this is the case we con¬ 
sider. Assuming the value of the ratio v'^l{gr) to be 
given from observations expression (I77D is basically a 
quadratic polynomial in Tb/T„ with parameters a and 
9. In figure [T] we show (1771) for 9 = 60° (upper panel) 
and 9 = 20° (lower panel), for different relevant values of 
a. To draw the figure we have used v^l{g r) = 0.06 (i.e., 
Vp = 5 km s“^ and r = 1.5 Mm), which is the value that 
led us in Sec. 12.11 Equation [T71 to conclude that there is 
no possible support for the plasma in the purely rotat¬ 
ing case. For magnetic support, the relevant stretches of 
the curves are those near the dashed horizontal line at 
ordinate = 1. For ease, we have marked the cut of each 
curve with that horizontal line with a large black dot and 
a thin dotted vertical line. The following guidelines have 
been used in the choice of values for a: in a solar tor¬ 
nado a is probably negative, i.e., the density decreases 
as one goes outward along the tornado field lines. Also, 
|a| should not be larger than of order unity, since the 
lengthscale of variation of p should be not too different 
from r itself. Finally, in the figure we restrict ourselves 
to concave-upward parabolas (i.e, a > —2), since those 
are the most favorable cases for magnetic support. 

We see that it is not difficult to find values of Tb/T^ 
that yield magnetic support of the plasma against grav¬ 
ity. In all cases shown, though, the basic variable |Tb/T„| 
must be above I: for 0 = 60°, for instance, the values 
marked in the figure range from —6.3 to —3.3, on the neg¬ 
ative side and from 4.3 to 7.3 on the positive side. Note 
that the actual ranges are [—14.2, —3.3] and [4.3, oo]. For 
9 = 20°, the absolute values are smaller. We expect Ty 
to be perhaps 1 or 2, reflecting the fact that the outward 
motion of the tornado is possibly as fast as the rotating 
motion or a little less so. Thus, for magnetic support one 
would need a substantial level of magnetic twist, possibly 
\Bp/Bp\ from a few to several units. If the poloidal flow 
becomes less important, then the amount of magnetic 
twist would increase to unrealistic values. 


3. CONCLUSIONS 

In this work we have investigated a possible mechanism 
to support dense plasma in a prominence barb against 
gravity. We have modeled a barb tornado as an axisym- 
metric structure with rotation about the symmetry axis, 
in which the magnetic field is vertical close to the tor¬ 
nado axis. Pressure gradients are ruled out as a support 
for the plasma because the small pressure scale height of 
prominence plasma implies that they could support ver¬ 
tical structures only a few hundred kilometers tall, much 
smaller than the height of observed barbs. 

In a barb tornado with a rotating helical field, extra 
magneto-centrifugal forces are present. We have found 
that the centrifugal force is much smaller than solar grav¬ 
ity for the barb-tornadoes observed so far. However, the 
poloidal magnetic force is a good candidate to support 
cool prominence plasma or even to inject such plasma 
into a prominence. For that, the structure must have sig- 
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nificant magnetic twist and/or poloidal flows not much 
smaller than the rotational velocities. Whether this is or 
otherwise the case is a question for future observations. 
However, more theoretical work and observational evi¬ 
dence are needed to elucidate the origin of poloidal flows 
or magnetic twist. 



Tb/T, 


Fig. 2.— Magnetic force along the field lines compared with the 
longitudinal component of gravity. For magnetic support, the value 
1 must be reached (marked by dots for the different curves) 
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